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1. Background 

The deconfining phase transition of SU(N) {N > 3) gauge theories in 3 + 1 dimensions is 
characterized by a low-temperature confined phase, where Z{N) symmetry is unbroken and quarks 
and gluons are bound, and a high-temperature deconfined phase where Z{N) symmetry is sponta- 
neously broken and quarks and gluons are free [1]. Simulations of pure gauge theories indicate that 
the transition between the confined and deconfined phases is first order for all A'^ > 3. The global 
Z{N) symmetry appears to always break completely in the deconfined phase, with no residual 
unbroken subgroup. 

The confined phase of pure gauge theories is in a region of low temperature that cannot be 
accessed perturbatively. It is therefore useful to generalize the system to restore the confined phase 
in a region of high temperatures. We were motivated in part by Davies et al. who generalized 
the mechanism of color confinement in a monopole gas to four-dimensional supersymmetric gauge 
theories on X 5' [2]. They showed that monopole contributions to the superpotential led to 
an effective action with a Z{N) symmetric minimum, corresponding to the confined phase, for 
all values of the 5' circumference (naively analogous to temperature). Therefore it is reasonable 
to expect that the addition of a term to the pure gauge theory action which mimics the effects of 
monopoles would make the confined phase to accessible at all j3. To this end we extended the 
Euclidean action of the pure SU (N) gauge theory with a Z{N) invariant term, the adjoint Polyakov 
loop: 

- I d^x Ha TrAPix) = -T dt j d^x Ha TrAP{x) . (1.1) 

Here P{x) is the Polyakov loop at the spatial location x, given by the path ordered exponential of 
the temporal component of the gauge field. 

A heuristic argument suggests that confinement is restored at high temperatures through vari- 
ation of /?A. Consider minimization of the effective potential 

Veff = £ VRrr^P - ThATrAP. (1.2) 

R 

I 1 2 

Because TkaP = \TrfP\ — 1, positive Ha favors maximization of TkaP, which implies \TrfP\ > 0. 
Thus Z{N) symmetry is broken which suggests this region is in the deconfined phase. Negative Ha 
favors minimization of Tr^P, implying TrpP = 0, which defines the confined phase. Therefore for 
sufficiently negative Ha, the confined phase may be restored above the normal Ha = deconfine- 
ment temperature. In the weak-coupling regime of high temperature we can calculate the effective 
potential, pressure, string tensions and 't Hooft loop surface tensions and examine their behavior 
in the restored confined phase resulting from the variation of Ha (see also [3]). 

2. SU{3) Simulation Results 

Our simulations were performed in SU (3) and SU{4) by adding an adjoint Polyakov loop term 
to the standard lattice action: 

S = Sw+Y,HATrAP{x) (2.1) 
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where Sw is the Wilson action. The naive relationship between the variable lattice parameter Ha, 
and the parameter used in our analytical calculations Ha, is Ha = hAO^, but there is an additional 
unknown multiplicative renormalization factor. 



Phase Diagram for SU(3) L=4x24x24x24 
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Our simulation results for SU{3) show 
that increasing positive Ha decreases the de- 
confinement temperature as expected, and 
for sufficiently negative Ha confinement is 
restored at high temperature. However, for 
negative Ha there is an unexpected new 
phase which breaks Z(3) symmetry in a pe- 
culiar way. 

Figure 1 shows the phase structure in 
the p - Ha parameter space of SU (3) defined 
in terms of (TrpP), where projection to the 
nearest Z(3) axis is understood. In the region 
of negative Ha there are 3 distinct phases: the 
deconfined phase with (TrpP) > 0, the con- 
fined phase with (TrpP) = 0, and the new 
"skewed" phase with (TrpP) < 0. 

As shown in Figure 1 , decreasing Ha at 
fixed /3 > 6, we encounter first the decon- 
fined phase, then the skewed phase, then the confined phase. To obtain the locations of the phase 
transitions we use the histograms of the fundamental Polyakov loop in combination with plots of 
the adjoint Polyakov loop susceptibility. Figure 2 shows SU{3) histograms of the fundamental 
Polyakov loop order parameter. 
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Figure 1: Phase diagram in SU{3) for an extended 
action 
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Figure 2: SU (3) Polyakov loop histograms 



(c) skewed, with tunneling between 
phases 



Figure 3 shows (TrpP) versus Ha- The presence of all three phases is clear. Figure 4 shows 
the adjoint Polyakov loop susceptibility. The obvious discontinuity in the order parameter shows 
that the transition between the deconfined and skewed phases is first-order in both graphs. The 
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0.4 




Figure 3: Projected TrpP Figure 4: Adjoint susceptibility Xm 

transition between the skewed phase and the confined phase is much weaker. It is likely to be first 
order as well, but this is not obvious due to the very small changes of the order parameter. 



3. SU{3) Theory 



To confirm our lattice results we studied the thermodynamics of the system using an effective 
potential adapted from the one-loop free energy density first evaluated by Gross, Pisarski, and Yaffe 
[4], and by N. Weiss [5]. Our modified expression is 



(3.1) 
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where the sum is over Matsubara frequencies &)„ = 2nnT . To locate the phases it is useful to write 
this as a function of the eigenvalues of the Polyakov loop: 
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In SU{3), it is sufficient to consider Ve/f for the Polyakov loop projected onto the real axis, P 
diag[\,exp{i(l)),exp{—i(l))]. Figure 5 shows that the effective potential finds all 3 phases. 



4. Comparison of St/ (3) theory to simulation 

We have calculated the values of (j) that minimize Ve/f in the three phases, then found the 
location of the phase transitions in terms of the dimensionless quantity Ha/T^ ■ The deconfined- 
skewed phase transition is located at Ha/T^ = —71^/48 ~ —0.206. The skewed-confined phase 
transition is at Ha/T^ = — 57r2/162 ~ —0.305. The ratio of these values is similar to that from 
simulations. 
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Figure 5: Phases from calculation ofVeff in 5t/(3) 
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Figure 6: Theoretical prediction for the pressure from 
Veff normalized by the black body value as a function 
of hji^. 



We also compared values for the pres- 
sure determined from the effective potential 
to the pressure determined from simulations. 
Figure 6 shows the theoretical pressure from 
the effective potential. From simulations the 
pressure is calculated along a path of con- 
stant j8 using 
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Nf / dHA{TrAP) (4.1) 



Comparing Ap/T'^ across the deconfined and 
skewed phases we find for theory Ap/T'^ = 
Ti'^/G ~ 1.64 across the deconfined phase, 
and Ap/T'^ = across the skewed phase. 
In simulations Ap/T^ = 1.64 it 0.03 across 
the deconfined phase and Ap/T'^ = —0.18 ± 
0.07 across the skewed phase. 



5. SU(4) Simulation 

The case of SU (4) is somewhat different. In simulations the new phase is partially-confining 
instead of skewed. Figure 7 shows the SU (4) histograms of the fundamental Polyakov loop. The 
new phase again occurs for negative Ha- We first encounter the deconfined phase, then the partially- 
confined phase. Tunneling is observed as we continue decreasing Ha in the partially confined 
phase. The fluctuations gradually reduce in size, but we are uncertain if there is a transition into 
the confined phase. 

In the new phase of SU{4), global Z(4) symmetry breaks spontaneously to Z(2), a partially- 
confined phase. The residual Z(2) symmetry ensures that (TrpP) = 0, but that {TvrP) ^ for 
representations that transform trivially under Z(2), so quarks are confined, but diquarks are not. 
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(P(.i)) in SU(4) ,3 = 11, fc.i = -0.13 
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Figure 7: SU{A) Polyakov loop histograms 



The Z(2) symmetry of the partially-confined phase is clear from the time history of variations of 
the real and imaginary parts of the Polyakov loop during a long run in which tunneling is observed, 
as shown in Figure 8. 
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Figure 8: Real and imaginary parts of SU{A) Polyakov loop versus Monte Carlo time 



6. SU{A) Theory 

For our analytical calculations in SU (4) we use again the one-loop effective potential to exam- 
ine the possible occurrence of four different phases: the confined phase with full Z(4) symmetry, 
the deconfined phase, a partially-confined Z(2) -invariant phase, and a skewed phase similar to that 
of SU{3). However, only the deconfined phase and the Z(2) phase are predicted by our simple 
theoretical model. A more complicated model with additional terms should reveal the confined 
phase [6]. 

We compared the phase structure predicted by the one-loop effective potential with our sim- 
ulation. Veff predicts a first-order transition between the deconfined and Z(2) -invariant phases at 
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Ha/T^ = -7i^/4S ~ -0.205617. This is in the same region of as in simulations. Across the 
deconfined phase, the theoretical value of A {^p/T^^ = 71^/3 ~ 3.289. In simulations A [p/T^) = 
2.21 ±0.07 

7. Discussion and Conclusions 

We have considerable evidence, from lattice simulation and from theory, for the existence 
of new phases of finite temperature gauge theories in SU{3) and SU{4) when a Z(A'^) -invariant, 
adjoint Polyakov loop term is added to the gauge action. In SU (3), confinement is restored at high 
temperatures, and the skewed phase was found. 

It is interesting to note that Wozar et al. [7], in their study of SU{2>) spin models, observed 
a number of interesting new phases. One of these, which they refer to as the anti-centre phase, 
appears similar to our skewed phase. The anticenter phase resulted from an action of the form: 

^eff = ^fSf +^\5S\^ (7.1) 

which includes a nearest neighbor coupling term in the 15 representation instead of an adjoint 
potential term. We believe that these phases are related. 

In the general case of SU {N) , there is good reason to expect a very rich phase structure may 
exist. For example, in SU (6), we can consider partial breaking of Z(6) to either Z(2) or Z(3). We 
have calculated the string tensions and 't Hooft loop surface tensions in the restored confined phase 
at high temperature [6]. These predictions can be checked in lattice simulations. 
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